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Abstract 



We analyze the gravitational wave signatures of a network of metastable cosmic strings. We con- 
sider the case of cosmic string instability to breakage, with no primordial population of monopoles. 
This scenario is well motivated from GUT and string theoretic models with an inflationary phase 
below the GUT/string scale. The network initially evolves according to a scaling solution, but with 
breakage events resulting from confined monopoles (beads) being pair produced and accelerated 
apart. We find these ultra- relativistic beads to be a potent source of gravitational waves bursts, 
detectable by Initial LIGO, Advanced LIGO, and LISA. Indeed, Advanced LIGO could observe 
bursts from strings with tensions as low as ~ 10^^^. In addition, we find that ultra- relativistic 
beads produce a scale-invariant stochastic background detectable by LIGO, LISA, and pulsar tim- 
ing experiments. The stochastic background is scale invariant up to Planckian frequencies. This 
phenomenology provides new constraints and signatures of cosmic strings that disappear long before 
the present day. 
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1 Introduction 

Cosmic strings can be produced in any phase transition where the vacuum group manifold is not 
simply connected, e.g., during the breaking of a U{1) symmetry. As topological defects, they appear 
to be completely stable at the level of the effective action where the U{1) symmetry is manifest. 
However, stability may break down in the ultraviolet. If, for example, the U{1) symmetry is 
embedded into a non-abelian gauge group at some higher energy scale, there must exist monopoles 
somewhere in the spectrum. The flux from these monopoles is confined under the U{1) symmetry 
breaking, and so hybrid defects can exist with strings ending on beads (we distinguish between 
monopoles and beads: the latter have their flux confined into strings). A second source of instability 
for global defects can arise from non-perturbative effects such as instantons, which generically 
lift the vacuum degeneracy to isolated points. The cosmic strings then become boundaries of 
domain walls. Both instabilities can lead to a rapid demise of the network, rendering cosmic 
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strings metastable at best. The question of cosmic string stability in the context of hybrid defects 
was originally investigated in Grand Unified Theories (GUT). The Langacker-Pi mechanism pTj, 
whereby monopoles created at the GUT phase transition have their flux confined into strings, was 
proposed to solve the monopole problem. It eliminated the monopoles by confining their flux into 
cosmic strings, which would in turn oscillate and quickly decay to radiation. In this paper, we 
will be interested in looking at the phenomenological signals of such a decaying network of cosmic 
strings. 

Cosmic strings formed in the early universe could result in a myriad of astrophysical phenom- 
ena. These include ultra-high energy cosmic rays j2], gamma ray bursts |3], radio wave bursts 
magnetogenesis |5], strong lensing P,[71[8], weak lensing |9], microlensing [101 [H], along with effects 
on the cosmic microwave background (CMB) polarization [12], the CMB spectrum at small angular 
scales [13l[Tl], and the cosmic 21 cm power spectrum [15] . 

One of the most exciting recent predictions is the observability of cosmic strings via their grav- 
itational wave signals. Cosmic string produce powerful bursts of gravitational radiation that could 
be detected by interferometric gravitational wave detectors such as LIGO, Virgo, and LISA [T6l fT7] . 
In addition, the stochastic gravitational wave background can be detected or constrained by vari- 
ous observations including Big Bang Nucleosynthesis (BBN), pulsar timing experiments [18], and 
interferometric gravitational wave detectors [TBI [TU] . 

An appealing source of cosmic strings from string-theoretic models such as brane inflation has 
led to a renewal of activity in recent years [20] . In such a scenario, the infiaton is represented by the 
relative location of a D-brane within the compactified extra dimensions. In some models, inflation 
ends with brane-antibrane annihilation, and the resulting open-string tachyon condensation leads 
to the formation of co-dimension two topological defects, as is typical in hybrid-inflationary models. 
In the generic case, a network of cosmic -F and/or -D strings is produced. The phenomenology of 
a network of stringy cosmic strings is richer than the standard abelian Higgs cosmic strings due to 
new parameters such as the spectrum of tensions [211 ES] , intercommutation probability [23j (see 
[21] for discussion of bursts from intercommutation) lensing [25], and as we shall discuss here, a 
decay rate per unit length. 

While the production of cosmic strings is essentially guaranteed following brane anti-brane 
annihilation (or more general hybrid inflationary models), the stability of such strings is model 
dependent. In the context of heterotic string theory, it was realized early on that the strings would 
be unstable due to the formation of domain walls stretching between various parts of the network 
|26j (the tension was also too high in these early models). The situation was reanalyzed in the 
context of Type IIB flux compactification first by Copeland, Myers and Polchinski |2I] and another 
source of instability, to breakage, was discussed. These two instabilities are complimentary, in the 
sense that typically exactly one will occur. This was further analyzed in [27] where it was shown 
that the generalized Green-Schwarz mechanism eliminates the domain wall instability leaving only 
the possibility of breakage on beads. The reason breakage is merely a possibility is because it is often 
not clear what the monopole is (or even if it exists) at the level of the low energy effective action. 
Multiple examples of stringy monopoles are known, ranging from Bl-ons [2S] to configurations of 
branes (for a recent example see [29]). In some cases, the strings carry a fraction of the minimum 
magnetic charge, and the network of strings, instead of having endpoints, has multiple strings ending 
on a given bead. (An important example is the "baryon" D3 branes described in ^30]). This leads 
to so-called cosmic necklaces, and they have interesting phenomenology of their own [HH E21 E31 El] • 

In this paper, we will be interested in the scenario where the strings are unstable to breakage due 
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to monopoles/beads. We will not rely on any specific string-theoretic models, and we will assume 
vanilla cosmic strings and beads; they are neutral under every long range force except gravity. 

Group theoretic considerations require that the symmetry breaking scheme must produce monopoles 
at a higher scale than the subsequent cosmic strings. The pattern of symmetry breaking is [35| (see 
also [36] for more discussion on cosmic strings in GUT models) 

G ^ H xU{l) ^ H , (1) 

where G is a semi-simple Lie group. Monopoles are created under the first symmetry breaking, 
and their flux is confined under the second, forming a network of beads and strings. A corollary of 
such a sequence is that the cosmic strings are at best metastable. Such a scenario has previously 
been considered EH] , and it was found that the beads rapidly annihilate due to the confining 
potential generated by the strings. The close proximity of beads means that such a network is 
unstable. However, in the special case where the monopole creating phase transition occurs before 
inflation and the string production phase after inflation, there will not be any primordial population 
of beadsj^ Such a network is merely metastable. Strings will then only decay via the Schwinger 
process (see |ES]), whose time scale is given by the decay rate per unit length 

r2 = ^exp (-vrmV/i) , (2) 

where n is the string tension and m is the bead mass. We parametrize the bead mass in units of the 
string tension by the parameter k = m? / ^. From the symmetry breaking pattern above k > 1. A 
derivation of Eq. (|2]) is given in Appendix |Aj For large values of n (for example, k > 86 if the tension 
is GyU ~ 10^^), the cosmic string lifetime is greater than the age of the universe, and the network 
is stable for all practical purposes. For small values of k, the network decays almost immediately 
upon formation. 

Martin and Vilenkin first studied the gravitational waves produced by these hybrid defects jl^ . 
They calculated the power radiated per frequency interval after solving the Nambu-Goto equations 
of the string/bead system. They assumed a network with a nonzero primordial monopole — >■ bead 
population and calculated the total power and stochastic background of the resulting gravitational 
waves. In this paper we extend their calculation to include an analytic solution of the waveform 
in different frequency intervals, and we consider the string motivated scenario in which the initial 
population of monopoles is absent due to inflation. 

We also calculate the phenomenology of highly focused gravitational wave bursts. In a celebrated 
paper [TB] , Damour and Vilenkin have shown that a network of cosmic strings can be observed via 
experiments like LIGO and LISA through the gravitational wave bursts produced by cusps and 
kinks on string loops. We find that intense bursts are produced by the acceleration of the ultra- 
relativistic monopoles, and we present the range of parameters where we can hope to detect these 
bursts. This calculation is of theoretical interest in its own right since a cosmic string ending on 
beads is the archetype of gravitational radiation emitted from a linearly accelerated mass, a problem 
famous for its simplicity and paradoxes, as we will review. 

Our main results are that for a moderate tuning of the monopole mass. Advanced LIGO will 
be able to detect strings whose tension exceeds ~ 10~^^. We also find that a stochastic signal 

brane inflation models, cosmic strings are naturally produced at the end of inflation, while monopoles are not 
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persists for networks which decay at any time between reheating and the present day. Because 
loop production is highly suppressed when the typical length of segment becomes sub-horizon, the 
gravitational wave phenomenology from an unstable network is complimentary to that of stable 
cosmic strings: an initial "stable string" epoch of loops is usurped by oscillating segments, followed 
eventually by the complete demise of the string network. 

The spectrum is very nearly independent of both frequency and k = m? / fi. The frequency 
cutoffs range from the Hubble scale to well beyond the Planck scale. This is an indication that 
massive {e.g., dilaton) radiation may contribute to the phenomenology |¥T| 

In ^ we review the issue of met ast ability, with particular emphasis on the recent models 
motivated by string theory. In ^ we calculate the gravitational waveform and power radiated by 
strings ending on massive beads. In ^ we track the cosmology of such a network and find the 
density and length distribution of cosmic strings. We then calculate the observational signals from 
such a network in ^ before concluding. We have relegated many of the more technical details to a 
series of appendices. 

2 Metastable Cosmic Superstrings 

Cosmic strings can decay by either becoming the boundary of a domain wall (the strings are 
confined) or by developing boundaries themselves through nucleation of bead endpoints (the strings 
are screened). Because the boundary of a boundary does not exist, these two types of instabilities are 
mutually exclusive. To see how they arise, let us start with a global cosmic string whose dynamics 
is well described outside the core by the Kalb-Ramond action 

S = [ + [ B2~ [ d^x ^H^^.H^^'^P , (3) 

where S is the string worldsheet, B2 is a 2-form with = dB2 the associated field strength, /i is the 
string tension and 77 its charge. We have suppressed the Einstein-Hilbert term. In 3+1 dimensions, 
the massless 2-form is dual to an axion: dB2 = -kdcj). The Kalb-Ramond action (Eq. [3]) therefore 
represents a cosmic string with a magnetically sourced long-range axionic force. There is a global 
[/(I) shift symmetry of the axion — + c, and since this string carries a global charge, local 
breakage is forbidden by topological considerations. Interestingly this fails to protect the strings, as 
domain walls will confine and rapidly eliminate them [l3l |2T] . This occurs when instantons generate 
a periodic potential for the axion of the form V{(j)) = A^(H-cos0) where A is some non-perturbative 
energy scale. At the minimum of the instanton potential, the global U{1) is broken to a discrete 
subgroup, and domain walls arise. As the axion winds around the global strings, it will necessarily 
have "to go over the bump" of the periodic potential; therefore the string is always going to be 
the boundary of a domain wall (see Fig. [T]). In general, the wall's tension can be high and will 
confine the strings together. Hence we expect a rapid demise of the string-domain wall network. 
The phenomenology of such an hybrid network may be interesting, and so we leave this for future 
work. 

We will be interested in the case of a local string which has no long range axionic force. The 
abelian Higgs model is the archetype of such a local cosmic string, but it is illuminating to minimally 
generalize the Kalb-Ramond action to include a gauge field. Consider the action 

S = -^ [ d^a + 7] [ B2- [ d^x {^^H,,pm^P + \F^,F^^-) + iB2 A F2 , (4) 
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Figure 1: As the axion winds around the string, the non-perturbative potential V{(p) must attain 
its maximum in some direction. A domain wall thus confines any string. Such walls will pull infinite 
strings towards each other and cause rapid demise of string loops. 

where F2 = dA is the U{1) field strength and the last term is the coupling between the axion and 
the gauge field parametrized by ^. This can be shown to be equivalent to the abelian Higgs model in 
the limit where the massive Higgs is integrated out of the theory. The B2 A F2 spontaneously breaks 
the U (1) gauge symmetry and the gauge field picks up a mass ^. This is just the Stiickelberg model. 
(In string theory, these types of terms can implement the generalized Green-Schwarz mechanism for 
anomaly cancellation). The massive gauge field has eaten the axion, and it now propagates three 
degrees of freedom. There is no long range axionic interaction, and the charge is therefore screened. 
Gauge invariance forbids an instanton potential for the axion, and it again appears as if the string 
is completely stable. 

This is no longer generically true when we consider the possible UV completion of the theory. 
Indeed, the U{1) could have been embedded in a non-abehan gauge group at some higher scales, 
in which case heavy monopoles will exist for the string to break on. In string theory, more exotic 
UV completions can also lead to monopoles. For example, by generalizing the kinetic term of the 
gauge field (which we took to be canonical in Eq. |4]) to the Born- Infield action (a non-linear gener- 
alization of Maxwell's theory), one can show the existence of monopole-like solutions representing 
the endpoint of a string on a brane [28j. That is, even if the gauge group remains U{1) at high 
energy, a cosmic (p, g)-string can still break on a brane in string theory and the endpoints will look 
like monopoles (or dyons) in four dimensions. The only way to forbid the existence of monopoles 
which could serve as a given string's endpoint is to use the Dirac quantization condition, whereby 
such a local string is properly classified as a stable Aharonov-Bohm string. 

Polchinski has conjectured that monopoles of minimum Dirac charge always exists in string 
theory jSj. This reduces the stability question to that of the existence of particles of fractional 
electric charge (in units of the abelian Higgs charge). If no fractional electric charges exist, than 
the Dirac monopole will carry the same charge as the string (fiux tube), and breakage can occur. 
If particles exist of electric charge 1/M, then bead nucleation will only permit g-strings to decay 
into q — M strings, and so only M coincident strings can break completely. One can thus look at 
stability in two equivalent ways. On the one hand, if the necessary magnetic charge for a string to 
terminate on does not exist, the string cannot break. To forbid the existence of such an object in 
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Figure 2: In the context of warped compactifications, different mass scales for tlie beads/monopoles 
are possible depending on the specific model. Heavy (bulk scale) beads (left) form when a distant 
throat is responsible for the instability. Light beads (k > 1) form (right) when branes exist near 
the tip of the inflationary throat. We assume inflation has already ended in these illustrations. 
The orange (dimmer) rectangle represents the far IR, where the strings reside. Black rectangles are 
D-Branes. 

a UV independent way, we simply require that the charge quantization satisfy Dirac's bound 

min ^1 ^min i (5) 

and then introduce particles of electric charge equal to a fraction of the abelian Higgs charge. 
Alternatively, one can show that by introducing particles of fractional electric charge, there will 
exist an observable Aharonov-Bohm phase upon transport of these particles around the string 
at arbitrary distances. This topological phase forbids breakage (assuming no additional massless 
fluxes exist) for the same reason that a global string cannot break, namely causality. Experiments 
far from the string should not be sensitive to a simultaneous event near the core. If no experiment 
can detect the enclosed string, breakage can occur, and so the lifetime of the strings becomes a 
phenomenological parameter. In summary, there are three different cases: 

• Global strings which are unstable to domain wall formation. 

• Local strings with no Aharonov-Bohm phase, these are not stable due to breakage via monopoles. 

• Local strings with an Aharonov-Bohm phase detectable at infinity. These strings are com- 
pletely stable. 

In this paper, we will exclusively discuss the second case where the strings break on monopoles/beads. 
The lifetime is going to be determined by the mass of these beads. Because this is strongly model 
dependent, we will simply illustrate a range of possibilities in the context of brane infiation in a 
warped throat [15] (see Fig. [2]). 
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In this scenario, the cosmic strings are formed at the bottom of a warped throat where the brane 
and anti-brane annihilate. If there are spectator branes near the tip after annihilation (which can 
be stabilized at various positions the string should be able to break directly on them and we 
expect the mass of the monopoles to be of the same order as the string tension (both set by the 
warped string scale), and therefore we expect k ~ although numerical factors can be quite 

important]^ If the other branes at the tip are slightly separated, then k can be greater than one, 
while if the branes are in different throats then we expect k ^ 1 and a very long-lived network 
(this case was emphasized in [21] )• More exotic monopoles configurations (from various wrapped 
branes) can also exist in these constructions |29]. We expect k > 1 for these monopoles where the 
numerical factor depend on the specific setup under consideration; for the rest of this paper we will 
take it to be a free parameter. We find that the decay of the network is observable for values of 
K < 80, which is not necessarily a fine-tuned mass range for the beads. 



3 Gravitational Wave Signal 

A straight cosmic string with two bead endpoints will oscillate back and forth once its length is 
subhorizon, producing gravitational wave bursts in the process. This is illustrated in Fig. A 
cosmic string with endpoints is the archetypical consistent system of piecewise uniform acceleratioE|^ 
The beads are being pulled by the string and accelerate with constant proper acceleration a = fi/m, 
in the straight string approximation. At the fly-by, the acceleration abruptly switches sign, and this 
discontinuity releases a burst of radiation of characteristic spectrum 1//^. Interestingly, this is not 
responsible for most of the radiated power. It is the ultra-relativistic oscillation which produces the 
most radiated power from the segments, as was first pointed out by Martin & Vilenkin [ID]. The 
spectrum from this dominant piece has the scale invariant 1// form. 



3.1 Radiation from a Uniformly Accelerated Mass 

Radiation from uniform acceleration is plagued with several conceptual issues [ITJ HHl SSI EOl ED [S2] ■ 
As in electrodynamics, the crux of the matter has to do with a radiation reaction which is pro- 
portional to X and therefore vanishes for an hyperbolic trajectory. This has led some physicists 
(such as Pauli) to incorrectly proclaim that there should be no radiation (in the case of electro- 
magnetism) from uniform acceleration, since it would violate energy conservation. The solution 
to Maxwell's equations for an isolated, uniformly accelerated charge is both simple and singular: 
the field strength on the light-front is proportional to a delta-function. It is known that a particle 
undergoing uniform acceleration over an arbitrarily long finite proper time does radiate. This fact 
does not undermine conservation of energy or the equivalence principle. We expect a similar result 
to hold for gravitational radiation. 

In the next section we will compute the spectrum of gravitational radiation from the piecewise 
uniformly accelerated beads following the work of [40]. This problem has already been studied 

^This was pointed out to us by J. Blanco-Pillado. 

^In E&M one can imagine an external force that provides the uniform acceleration. In GR, this "external force" 
will perturb the metric as well and should be included. The solution for a uniformly accelerated mass in GR leads 
to a class of metric called C-metric [53] (see also [54]). This class of metric always include conical deficit angles, i.e., 
cosmic strings. Therefore, a cosmic string with endpoints is the self consistent way of getting uniform acceleration 
in general relativity. 
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Figure 3: An open cosmic string with bead endpoints. The beads are uniformly accelerated until 
the fly-by occurs at some minimum impact parameter b, at which point the acceleration abruptly 
change sign. The presence of kinks (not shown) on the string will also lead to abrupt change in the 
acceleration. 

extensively in the context of E&M (see for example [SS])- Before going into the detailed calculation, 
we can summarize all the key features. The waveform is a strong 1/f burst with the following 
properties: 

• The power radiated is approximately independent of the length of the segment 

• The power radiated is approximately independent of the mass of the bead 

• The radiation is scale invariant for 1 < // < 7q 

The scale invariance of the burst of radiation is a rather unusual feature. This should be thought 
of as a consequence of the ultra-relativistic kinematics. (A similar behavior was found for the 
ultra- relativistic rotating rod by Martin & Vilenkin.) Naively, one would expect the discontinuous 
acceleration to produce a strong 1//^ burst. Instead, the even stronger 1/f waveform is dominant 
whenever 70 ^ 1. 

3.2 Gravitational Waveforms 

For a periodic source with period T, the solution to the linearized Einstein equations in the wave 
zone has the form 

00 

V(x,t)= e!rj(x,c.„)e-''"-% (6) 

n=— 00 

where 

4-C 1 

e(:)(x,^„) = — (T,.(A;„) - -r/,.T,^(A:„)) , (7) 
is the polarisation tensor, and 

T^'^iK) = fj^ dt j rf3xT^^(x)e^^"-^ , (8) 



8 



is the Fourier transform of the stress-energy tensor. The wave vector is given by kl^ = c<j„(l,f^), 
with Un = 27rn/T, and f2 is a unit vector pointing from the source to the point of observation. If the 
bead-string system is confined to the 2;-axis then the stress-energy tensor only has three non-zero 
components: T°°, T^^ and T^^. The conservation equations V^T^"^ = in Fourier space is 

UnT^"" - k'T" = , (9) 

which leaves only one independent component. 

We need the stress-energy tensor of the idealized setup of a straight string connecting two beads. 
This was first calculated by Martin & Vilenkin [50] directly from the Nambu-Goto equations of 
motion. The stress energy tensor (simplest choice is T^^) is given in the time-domain by 

T°3(t,x) = m{^oVo-a\t\)[5{x-Xi{t)z)-5{x + xi{t)z)] , (10) 

where 

;sgn{t) 



x,{t) = (-1)^:^^(70 - + ilovo - a\m (11) 
a 

is the position of the two beads (by symmetry Xi{t) = —X2{t)) and 70 and fo are the maximum 
gamma factor and speed of the beads, respectively. In this case, the period of motion is T = 
where a = fi/m is the proper acceleration of the beads and the frequency is 

nira , , 

^n = , 12 

The maximum length of the system, when the beads are at rest, is / = and so for ultra- 

relativistic beads we have, to good approximation, that T ~ /. The Fourier transform of Eq. (10) 
is, 



T°^(u;„,k) = mjoVoIn{u) , (13) 



with 



In{u) = f idiicosin-Kil + u^e + l/(7o^o)2)) 

Jo ^0 



and u = k^/uJn- We may take the unit vector that defines the direction of travel of the wave to be 

cos sin 6 

n = I sin0sin^ | (15) 
cos 9 

so that u = cos 6. The Fourier transform of the stress-energy tensor is thus 



/ u 1 \ 





\ 1 1/n / 



(16) 
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and its trace is ^ 

Txikn) = m-foVo{u • (17) 

u 

The polarisation tensor reads 

e}.AX,^n) = ln{u) 

( u + l/u 1 \ 

u-l/u 

u-l/u 

\ 1 u + l/u J 



As shown in Appendix [C} the gravitational wave is linearly polarized and the plus polarization can 
be written as 

e+(x,t^) = In{u) ■ (19) 

r u 

The waveform in the frequency domain is simply obtained by multiplying by the period h{f, r) = 
Te+(x, cu). As was shown in |10], the large frequency {n ^ 7q) behavior of the waveform is 1//^. 
In Appendix |D| we re-derive this result using a continuous approximation that shows clearly that 
this tail of the waveform is coming from the discontinuity. 

Almost all of the energy is radiated in the frequency range 1 < n < 7q. A good approximation 
of the waveform in this regime can be obtained by dropping the second term in the square root of 



Eq. ( 14 ) and performing the integral 



^ _ 16G'm7ot;o/sin^(n7r(l -u)/2) 
r n^7r^(M"^ — 1) 

Changing variable to = r27r(l — u)/2 gives 

4Gm-foVol sin^^ 
h{f,r) = — — . 21 

The burst is highly focused around ^ ~ which corresponds to small ^. Together with the fact that 
n is large, we can drop the second term in denominator. All the angular dependence is therefore 
describe by ^^y-^ which is zero at ^ = [0, vr, ■ ■ ■ ] and has its first maximum at ^ = \/2. The beaming 
angle can be well approximated by the position of the second zero at ,^ = vr 

0,^^, (22) 

'n 



This can be written in term of the frequency using n = f{l + z)l. We will consistently refers to / 
as the observed frequency, while / is the physical (maximal) length at time of burst. The waveform 
at its maximum value is 

M/,r)^C^l (23) 
r J 
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with C = 2i/2sm!_y^^ This is vahd for 70 ^ 1, fo ~ 1. Note that there is no redshift factor (1 + z) 
in this expression. The waveform swiftly goes to zero for / < 1/T ~ 1// and the approximation we 
have made breaks down for n ^ Jq. Hence most of the power is radiated in the frequency range 

< / < • (24) 



(1 + ^)/ ' {l + z)l 

The 1// dependence of the waveform signals a very strong burstj^ and leads to a scale invariant 
spectrum with constant power emitted per logarithmic frequency interval, as we will now show. 

3.3 Energy and Power 

The power radiated per frequency interval is 110] 

Pn = 2n^7i'^Gfi^ [ du {l/u - uf \In{u)\^ , (25) 
Jo 



where /„(«) is given by Eqn.(14). Again, there are two interesting frequency intervals, the 'scale 
invariant' range 1 < n -C 7o, and the 'convergent range' 1 ^ 79 ^ n. We can analytically 
determine the power radiated in each mode in the limit 70 00: 

= 1^ (37^ _ 4Ci(27rn) + 01(4?™) + \og{An^n^) + 4n7rSi(27rn) - 2n7rSi(47m)) , (26) 
where Ci and Si are the cosine integral and sine integral functions, respectively, and •je is the Euler 



constant. We can find the beaming angle by examining the integrand in Eq. (25). In the 'scale 
invariant' range, the radiation is beamed in a cone of apex angle 9n ~ >/lJn, and the power in the 
^th j2iode is approximated by 

P„ « ^ . (27) 

n 

The 'scale invariant' part of the radiation is thus constant per logarithmic frequency interval, and 
it dominates the total power radiated. 



1 

P ^ AGfi^ XI " ~ ^ ln(7o )G'/i2 . (28) 



n 

n=l 



There is a small correction to this formula coming from the part at high frequency but it can 
be safely neglected for large enough 70. It is instructive to re-derive the waveform from energy 
considerations. For a single burst, the meaningful quantity is the energy 

dE ^ ^^h^ 1/(1 + z)l< f < 7oV(l + z)l , (29) 

^We can roughly define a burst as a wave with a long tail in the frequency domain. The longer the tail in frequency 
domain the sharper the emission was in the time domain. For comparison, bursts from cusps and kinks go like f~'^^^ 
and f~^/^ respectively PJj. The bursts we consider here are much stronger. 
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which is beamed in a cone of angle ^/ ~ a/1/(1 + z)fl. The gravitational wave energy from an 
isotropic burst is given by 

dE = ATTr^Mlfh^{f)df . (30) 
Thus if all the radiation is contained in a cone of solid angle ^j, the waveform must satisfy 

dE = r^e)Mlph\f)df , (31) 



and so comparing with Eq.(29) yields the waveform (assuming the burst is aimed directly at us) 



h(f)~^!f. (32) 

r J 

For observation within a narrow frequency band, the beaming angle appears narrowed by redshift. 
This is because an observation of a frequency / today corresponds to a higher emitted frequency, 
and thus a narrower cone of radiation. Since the energy in a gravitational burst must scale like 
1/a, all of the redshift dependence is in the observed beaming angle, and the waveform remains 
unchanged. This means that sources at very large redshift are just as bright but are less often 
directed at us. 

3.4 Beyond the Straight String Approximation 

In the idealized straight string approximation, at the time of fly-by all the energy of the string/beads 
system is located in a region smaller than the corresponding Schwarzschild radius, and a black hole 
should form. This will not occur for segments with some minimum impact parameter, due to 
angular momentum. Given that the Schwarzschild radius is r<j = , and that the energy of our 
system is of the order /i/, we will form a black whenever the impact parameter h satisfies 

b<-^^G^il. (33) 
ml 

In this paper we work under the assumption that for most strings, the impact parameter is always 
larger than this bound. This is legitimate, since for a network of scaling cosmic strings the segments 
will contain a number of kinks. Kink-bead interactions are very similar to fiy-bys, since in both 
cases the bead's acceleration is piecewise constant with sharp discontinuities. Such strings will not 
form black holes since the bead impact parameter is of order the length of the segment. If the 
average number of kinks per segment falls significantly below one, black hole formation must be 
considered. This may yield interesting phenomenology in its own right. 

One may wonder whether the presence of kinks will alter the observable signal from such a 
network and here we argue that it does not. We will assume that the kink angle is quite small, in 
which case it may be thought of as a perturbation on a straight string. 

Under this assumption, we expect a string with k weak kinks to produce the same total power, 
but divided into k separate bursts. Starting fully extended and at rest, the beads will accelerate 
toward the segment center of mass, and encounter k kinks on the way. The first kink will cause 
a burst in the low frequency range, since the bead has not had a chance to reach its top speed. 
Successive bursts will be higher and higher in frequency, since the bead is moving faster and faster, 
until the maximum speed is reached. 
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Since the kinks have httle effect on the bead velocity we can simply divide up the (ffat) energy 
spectrum of the straight string scenario into k distinct domains. This means that the each of 
the k bursts represents a different interval in frequency space, but their amplitude is identical to 
the simple case. In essence, the kinks act like a prism, dividing the burst into coequal frequency 
domains. For experiments sensitive to either a narrow frequency band or total power, there can be 
no distinction between a single co-linear bursts, and prismatically scattered bursts. 



3.5 The Rate of Bursts 

Let the number density of bead-antibead pairs connected by a string of length between / and / + dl 
at time t be 

n{l,t)dl. (34) 

By length here we mean the fully stretched out length of the system when the beads are at rest. 
Each bead-antibead pair oscillates in a time T ~ / and the number of bursts per unit spacetime 
volume we observe at time t is 

u{l,t)dl jn{l,t)dL (35) 

We have restricted ourselves to the case when the beads are relativistic and the radiation is highly 
beamed. We do not observe all bursts, rather we only observe a fraction 

A(^, /, ^) - i(I^^0(7o - l)e ((1 + z)fl - 1) (70^ - (1 + z)fl)) , (36) 

with 7o ~ fil/m. The first term on the right hand side is the beaming fraction for an angle 
9f ~ ((1 + z)fl)~^^'^, the first Heaviside B-function (second term) ensures that we are in the 
relativistic regime 70 ^ 1, and the second and third Heaviside G-functions ensure that we are in 



the regime given by Eq. (24) where the waveform Eq. (23) is valid. Although in Eq. (36) we have 
included the upper bound on the frequency for clarity, it turns out to be irrelevant. The condition 
7g > (1 + 2;)// can be turned into a bound on the length of the segment of string connecting the 
beads, / > (1 + z)f{Gfi)~^nlp, where Ip is the Planck length. For typical values of Gfi and k even 
for the largest possible redshifts the lower bound on / is within a few orders of magnitude of the 
Planck length. 

Time and distances can be written as cosmology dependent functions of the redshift t{z) = 
HQ^Lpt{z), r{z) = HQ^Lpr{z) and dV{z) = HQ'^(pv{z)dz. These functions depend on cosmological 
parameters and need to be computed numerically as described in Appendix A of [17], where a 
vanilla A-CDM model was used with recently measured parameters that can be found in ■ One 
can obtain reasonable (good to about 20%) analytic approximations for the cosmological functions 
using variations of the expressions introduced in [T6|. In particular the following expressions, 

= {i+zr'/\i+z/z,,r'/', (37) 

= l + z/3.5042' 

12^2(1 + ;,/i.6)-W2(i + z/z,,)-'/', (39) 



seem to work well (with Zeq = 5440). We note that for consistency, the value of Hq that should be 
used with these functions is 73 km/s/Mpc [56]. This allows us to write the rate of bursts we expect 
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to see from a volume dV{z) (the proper volume in the redshift interval dz), from bead-antibead 
pairs with lengths in the interval dl as, 

(7/? 

(l + ^)~V(/,z)A(Z,/,z), 



dV{z)dl 
d¥} 

— = H^'{l + zr'M^Ml,z)A{lJ,z). (40) 

The first {1 + z) factor comes from the relationship between the observed burst rate and the cosmic 
time. To compute the rate of bursts an estimate of the minimum detectable amplitude A^i^ needs 



to be performed (see section 5.1). A matched filter search returns the amplitude A of a template, 

h{f) = At{f) (41) 
where in this case t{f) = f~^, and the dimensionless amphtude is given by 

A = , (42) 

where again C = v 8 sin {V2)/n. We can write the length as a function of the amplitude, 

and change variables to write the rate as a function of the amplitude. In the end we perform the 
integrals to compute the rate above the minimum detectable amplitude 

and see how that rate compares with one event per year (say). In the following sections, we will 
compute n{l,t) and Amin and the results for the rate is shown for LIGO, Advanced LIGO, and 
LISA in Fig.Q. 

3.6 Stochastic Background 

A stochastic background of gravitational waves is produced by the incoherent superposition of all 
bursts ever produced by the network. It is typically expressed as the ratio of the energy density 
of gravitational waves to the critical density today. To derive an expression for the stochastic 
background we begin with the expression for the plane wave expansion for the metric perturbation 

/oo P 
df dnhM,^)e''^^^'~^-'^/''^ eim. (45) 
^ -oo J 

Here A = +, x labels polarisations, Q is the direction of propagation of the gravitational wave and 
ef^(r2) is the polarisation tensor. If the stochastic background is isotropic, stationary, unpolarised, 
and uncorrelated the ensemble average of the Fourier amplitudes /;,^(/, fi) satisfies [57j 

{h\if,n)hA'if\n')) = 6\n,n')6AA'Sif - n H{f). (46) 
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This expression defines H{f). The relation between H{f) and Qgyj{f) is 



H{f) = 1^ \fr figwd/i) . 



(47) 



where Hq is the Hubble parameter which we take to be Hq = 73 km/s/Mpc [56]. In the following 
we describe how to compute H{f) using Eq. (46). We start with an expression for the metric 
perturbation as a function of the frequency evaluated at a particular point in space (we take x = 
for convenience), 



(4J 



From Eq. (48) it is clear that if we multiply both sides of Eq. (46) by e^,(fi'), integrate over the 
sphere (J^a dfl') and sum over polarisations {J2a')^ obtain, 

{h*M\ Q)h^\f, 0)) = etmif - n H{f). (49) 

Similarly, multiplying both sides by ef^(fi), integrating over the sphere {jg2 dfl) and summing over 
polarisations {J2a)^ obtain, 



iKbif, 0)h''\f', 0)) = 16n6if - n H{f). 



(50) 



Here we have used the fact that ^ab^^A ~ 4- To remove the delta function we borrow a trick from 
|16j . In the limit f ^ f the delta function takes the formally infinite time interval 6{f — /') — > T, 
so that 

(|M/,0)r) = 167rriJ(/). (51) 

The factor of T will cancel eventually. We have seen that the gravitational wave signal is linearly 
polarised which means that for a wave traveling along the z direction we can write the metric 
perturbation as 

cos 2iIj sin 2iIj 
hab{f,0) = h{f) \ sin 2^ -cos 2^ |. (52) 





If we average over the polarisation angle ip we get (/i*{,(/, 0)h°'''{f, 0)) = 2{\h{f)\'^), and hence, 

1 



Thus we can write. 



H{f) 



f^gwd/l) 



8nT 



mm. 



(53) 



(54) 



The strain (Eq. 23 ) and the rate of bursts (Eq. 40 ) enter the ensemble averaged strain through 

dR 



{\him=T / dz / dih\f,zj) 



dzdl 

Thus, the stochastic background produced by the incoherent superposition of bursts is 



figw(/) 



47r2 

ml 



f dz dlh\f,z,l) 



dR 
dzdl ' 



(55) 



(56) 
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where the integrals are over all redshifts z and over all lengths / connecting the bead pairs. 

As pointed out by Damour and Vilenkin [16] , if we are interested in the confusion noise produced 
by the network of beads and strings we must avoid the biasing of f2gw(/) that could result from 
including large amplitude rare events. We can accomplish this by excluding events that occur less 
frequently than the typical duration of bursts. The duration of a burst as seen in an instrument 
with peak sensitivity at frequency / is 0{f^^). Only bursts that overlap, i.e., small amplitude 
frequent bursts contribute to the Gaussian background we are trying to estimate. We can follow 
the procedure developed in [19j . We begin by expressing the rate as a function of the strain and the 



redshift rather than the length and the redshift. We can do this change of variables using Eq. (145) 
to write 

dR I dR 

= (57) 

dhdz h dldz 

We then integrate the rate over all redshifts and find the strain for which the rate at and above 
that strain becomes smaller than the peak experiment frequency /. Namely we find the value of 
such that, 

poo poo njD 

R{> h,) = dh dz^j^, (58) 

Jh, Jo 



dhdz' 



and evaluate 



Equation (59) gives the so-called confusion noise produced by the network, a Gaussian background 
that needs to be included along with other sources of noise to estimate the total noise present in 
our detectors. If the confusion noise is sufficiently large it can have adverse effects on our ability to 
detect individual bursts. 

For stochastic background searches that use cross correlations between detectors, such as those 
performed by LIGO, it is not necessary to remove the large amplitude events in our estimate 
of figw(/), and Eq. (56) may still be used. Large infrequent bursts will contribute to the cross- 
correlation in our estimate of i7g„(/), even though cross-correlations are not the optimal procedures 
for finding such bursts. These large infrequent bursts are optimally detected using other techniques: 
Matched-filters if they are sufficiently large amplitude (as described above), or methods optimized 
for non-Gaussian stochastic background detection |58] . 



4 Cosmological Evolution 

The cosmic string network forms at the end of brane infiation (or similar hybrid infiationary scenario) 
without any beads present. This infinite string network begins to decay immediately, resulting 
almost exclusively in super-horizon length segments. As long as the typical segment length is 
super-horizon, the evolution will be identical to that of an infinite string network, and scaling is 
readily expected despite the presence of a small number of sub-horizon (oscillating) segments. 

At a time t^, the typical segment length becomes sub-horizon. We will call this the beginning 
of the short string epoch. Much like matter-era cosmic strings, these segments are smoother and 
straighter than the segments during the long string epoch. The strings during this epoch will 
cease to obey a scaling solution. This is because the segments will have a drastically reduced 
intercommutation probability due to low peculiar velocities. The segment energy density will thus 
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increase relative to the scaling solution during the short string epoch until either it begins to 
lose significant energy to gravitational radiation, or the universe becomes matter dominated. The 
violation of scaling during the short string epoch can thus increase the energy density relative to 
long cosmic strings by a factor of at most (G/i)"^/^, which is significant in our parameter space. The 
network will continue to fragment into shorter segments, causing the energy loss to gravitational 
radiation to eliminate most of the network around the time t = t^^. The primary phenomenological 
quantity that we must determine is the number density of string segments of a given energy. 



4.1 Number Density of Segments: n{l,t)dl 

Let n{l, t)dl be the number density of bead-antibead pairs connected by a string of length between I 
and I + dl at time t. As before, we define the length I of a segment by its maximum possible length: 
I = E/^. Then the energy density in cosmic string is given by 

n{l,t)nldl = Pes- (60) 

Notice we do not need to distinguish between "infinite strings" and finite segments: the former are 
more accurately characterized as having a super-horizon length. We do not track loops explicitly. 
Instead, we note that loops are not produced by segments which have begun to oscillate. Before we 
perform a careful derivation of n{l,t)dl, it is enlightening to first consider a simple estimate that 
can be done in the scaling regime. We can treat the cosmic string network as a single worldsheet 
where decay events are uncorrelated points on the sheet. This means a metastable cosmic string 
of decay rate per unit length r2 will have decay events obeying a Poissionian distribution. The 
probability of a worldsheet proper area A containing k nucleation events is given by 

Thus at a time t, the probability for a given string to be of length / is given by the probability 
that exactly zero decay events have occurred in the worldsheet area It, i.e., in the entire past of 
the segment, times the probability that one or more decay events did occur in an area tdl. (We are 



asking for a segment of length at least / and at most / + dl.) Using Eq. (61) the probability of a 
segment having such a length is given by 

P{1, t)dl = T2texp{-T2lt)dl . (62) 



We can relate P{l,t)dl to the number density n{l,t)dl via Eq. (60) 



nU, t)dl = pa, t)dl ^''^^l , „ (63) 
^ ' ^ ^ ' ^ P{V,t)pl'dV ^ ' 

where pcs{t) ~ p/t^ is fixed by the scaling solution. Thus we have 

nil, t)dl = Tl exp (-/tEs) dl (64) 

This simple argument is expected to be valid in the scaling regime when most of the strings are 
super-Hubble in length. The number density of segments of finite length goes to zero both when 
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r2 — > (the network is stable) and when r2 ^ oo (the network decays immediately). For a fixed 
length and time, n{l,t) has a maximum at r2 = 1/lt. Therefore the time (t^:) at which a significant 
fraction of the segments in the network are subhorizon / ~ t is 

(Z) ~ ^ t = U ■■= T-'^^ (65) 

To correctly account for the deviation from scaling and the gravitational back-reaction, we need a 
more thorough analysis. Because n{l,t) is a length density it satisfies the continuity equation 

dn{l,t) d /• ,j .\ d 



dt dl 



(jn{l,t)^ -3\{l,t) + g , (66) 



where gdldt is the number density of pairs produced with lengths between / and I + dl, at a time 
between t and t + dt. The function I incorporates the effects of the geometry on individual segments, 
i.e., Hubble stretching and drag. We can characterize explicit terms within g as being one of three 
types: loop producing, segment intercommutation, and segment breaking. Loop production provides 
a stable solution for pcs{t) = J n{l, t)ldl, the scaling solution. Two-segment intercommutation will 
not be explicitly important and breaking will eventually destroy the scaling solution. We write 

g = 5'loop + gic + 5'break • (67) 

In Appendix |B} we discuss each term in more detail. Here we simply argue that intercommutation 
and loop production ensure scaling in the absence of breaking]^ To make the problem tractable, 
we will coarse grain over the scale of loop production and treat it as a continuous shortening of 
segments. The result is an / which now includes the effects of both the geometry (expansion) and 
loop production, 

/ = iH + Loj. (68) 
21 

= 3HI — — (scaling solution) (69) 



It is readily seen that the Boltzmann equation (Eq. (66)) with g = and the above I has a scaling 
solution with pcs ~ p/t"^ in both the matter and radiation eras. We are now ready to include 
breakage with 

^breakt^^ = (^^ n{l' ,t)dl' -ln{l,t)^dl (70) 

The breakage formula can be motivated as follows. The breakage rate of strings of length / is r2l, 
which explains the last term. The first term can be understood by considering the process by which 
n{l, t)dl increases. This is entirely due to breaking longer strings, those of length V > I. The rate at 
which longer strings break to produce those of length between / and l + dl is given by the number of 
longer strings present (hence the integral) , and the measure of string where a break yields a shorter 
string of length between / and / + dl, i.e., 2dl. The factor of two is because the breaking point 



^The main effect of intercommutation is to create a small scale structure on the strings that leads to the creation 
of loops. Loop production is the dominant channel for energy losses in the network. It is also sufficient for scaling - 
one can turn off gravity and the "infinite" network still scales. 
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can be closer to either of the two endpoints of the longer string. Because string breakage should 
not affect the string energy density, we can easily check that ^ filn{l,t)dl is independent of 
which it is. Finally, the last ingredient we will need to solve the Bolztmann equation is the energy 
loss in gravitational waves as sub-horizon strings oscillate. We already calculated the total power 
radiated by an oscillating segment [iO]. 

P^41n(7o')G/x^ (71) 

and the energy comes out of the length of string connecting the pair (we think of the beads coming 
to rest, and with each oscillation the string gets a little shorter because of GW emission). Thus the 
length of string (when the beads are at rest) decreases as 

4„^-41n(7o')G/i 1<H-' (72) 

We will parametrize this as /g„ = —TGfi with typical values F ~ 50 due to large 7 factor of the 
beads. The phase transition at t = t^, leads to a network that consists of subhorizon segments, 
which are straightened by expansion, but not stretched. These will oscillate and fragment until 
their GW radiation becomes dominant, at a time of order t*^,, after which the network disappears. 
We can estimate t^^ by comparing the power lost to scaling versus the power lost to gravitational 
radiation: 

Jti^) ^ ^G^,' j^nil^m (73) 
^ = tJ^/TGjl = l/VrsFG/i (74) 

Notice t^:^: does not depend on expansion, unlike t^,. We should therefore think of 1^,^, as the 
network lifetime irrespective of the scaling solutioEpl 



4.2 Solutions 

During the long string epoch, the network will satisfy the scaling Boltzmann equation (now an 



integro-differential equation given gbreak in Eq- 70 ) 
dn{l,t) 



dt 



d_ 
dl 



3HI 



21 
1 



n{l,t) 



- 3Hn{l, t) + 5(break 



whose solution is 



n{l,t) = 4F^exp(-2F2/t) , 



(75) 



(76) 



where we have assumed radiation dominance for now. This solution is, up to a factor of two in 



the exponential, exactly what we obtained in our simplistic derivation of Eq. (64). After the long 



string epoch, the strings will begin to oscillate, which will drastically alter their appearance and 
behavior. In particular, we can now ignore Hubble stretching and loop production, and therefore 
the only contribution to / is from energy loss to gravitational radiation. The Boltzmann equation 
then reads 

dn{l,t) _^p^^^^;^^)_3^^(;^^^^^^^^^^^ ^77) 



dt 



dl 



^We thank K. Olum for pointing this out. 
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which has the solution 

n{l,t) = r^exp(-r2/t)W-exp(-|rG/ir2t^) . (78) 

V ^* 

where one recognizes the previously defined t^,^, = 1/ ^y^Gfl^2 in the exponential. The complete evo- 
lution of the segments (including the cases where the network decay in the matter era) is extremely 
well approximated by 



n{l,t) = "^^^ "^^^i^*' ^eg}} ^-r2Ut+min|f,t,,fea})^-|tVtL ^ (^79) 

where med{} is the median. The exact solution is given in Appendix [b] 

5 Results 

In this section we examine the parameter space of the theory to find detectable bursts and stochastic 
background signals. 

5.1 The Minimum Detectable Burst Amplitude 

To produce an estimate of the minimum detectable amplitude we follow the analysis in [T7| for 
cosmic string cusps closely. As we have seen, the gravitational waveforms in the frequency domain 
are 

h{f)=At{f), (80) 
with t(/) = /~^. We will use the conventional detector-noise-weighted inner product 



/■oo 

{x\y) = 43? / df 
Jo 



where Sh{f) is the single-sided spectral density defined by {n{f)n*{f')) = ^6{f — f')Sh{f), where 
n{f) is the Fourier transform of the detector noise. The template t{f) can be normalised using the 
inner product o"^ = (t\t) so that t = t/a, and {t\t) = 1. For an instrumental output s(t), the signal 
to noise ratio (SNR) is defined as, 

p^is\i). (82) 

In general, the instrument output is a burst h{t) plus some noise n{t), s{t) = h{t) + n{t). When the 
signal is absent, h{t) = 0, it is easy to show the SNR is Gaussian distributed with zero mean and 
unit variance. When a signal is present, the average SNR is 

(p) = mi)) + {{n\i)) = {Aai\i) = M (83) 

and the fluctuations also have unit variance. This means that for an SNR threshold pth, on average, 
only events with amplitude 

Ah > — , (84) 
a 
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Figure 4: Burst detection rate for LIGO (left), Advanced LIGO (center), and LISA (right). The sohd 
curve represents the BBN constraint on string tension from loops i.e., long-lived strings, whereas 
the dashed curve represents the BBN constraint on string tension from bead induced bursts. 



will be detected. The form of the Initial LIGO noise curve is well approximated by [60j, 




L09 X 10 



4, /100Hz\^ 
+ 1.44 X 10"^^ ( ) 



+ 1.28 X 10"*° 1 + 



^5) 



The first term arises from seismic effects, the second from thermal noise in the optics, and the third 
from photon shot noise. A reasonable operating point for the pipeline involves an SNR threshold 
of Pth = 4. The value of the amplitude is related to p, the SNR, and a, the template normalisation. 



= m / df 



ShU) 



df 



f- 



Skif) 



(86) 



Because of Eq. (84) an event with an SNR close to the threshold, p* ~ 4 will on average have 
an amplitude A^, ^ 6 x 10~^^. Events on average will not be optimally oriented so we need to 
increase the amplitude by a/5 to account for the averaged antenna beam pattern of the instrument. 
So we expect A^^^ ^ 10~^^ for the Initial LIGO noise curve. For Advanced LIGO we expect to 
do about an order of magnitude better so we ^ 10 For LISA, our estimate of the 

minimum detectable amplitude is Amin = 4 x 10"^^. This estimate was made using a LISA noise 
curve that includes confusion noise from galactic white dwarf binaries, assuming that some of the 
binary signals can be fitted out[ 
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5.2 Burst Detection by LIGO, Advanced LIGO, and LISA 

We want to calculate the rate of bursts 

poo poo 77-) 

R^A^,^ = dA dz— (87) 



^min ^ 



dzdA 



with 



dR (pv{z)jy{A,z)A{A,z) m {1 + z)ipr{z) 



dzdA Hl{l + z) 2i2\J BGfiHoA 



which can be obtained from Eq. (40) by doing a change of variable from / to A using Eq. (43). 
The integral over amplitude goes to infinity to good approximation. More precisely, there is an 
amplitude (length) cutoff set by the size of the horizon, since only strings with length I < 1/H will 
oscillate and produce bursts. 

With some approximations, it is possible to obtain an analytic expression for the observed burst 
rate. First we will work at large redshift z ^ Zeq to simplify all of our cosmological functions. The 



number density of segments in the radiation era is given by Eq. (78) and we can make the further 
approximation that the network decay suddenly at time t = t** and set the density to zero beyond 
that point. Therefore we divide the density into the scaling regime t < t^, the non-scaling (short 
string) regime < t < t** and the post string regime t^^: < t. 

(I t <U 

n(/,t) =4r2e-2r2'* J i_^v,u^^^ < ^ < (gg) 

[ u^<t 

With this density, we get 

-^OC-^<J (M!!;!<!g-9W4e-aA/(2.2) z (90) 

dzdA A^f I //5/2^8 e e > z > z„ \ j 

z z^-^ 

a oc . (91) 

The integrals over A, z can be done exactly, and the results can be written in terms of incomplete 
gamma functions. We found that for most parameter values, the detectable signal is dominated by 
the radiation emitted at z** which can be before or after Zeq- We therefore had to resort to numerical 
integration to incorporate the signal coming from low redshift. For all numerical work, we assumed 
the scaling epoch describes an average of 40 long strings per Hubble volume. The results are shown 
for LIGO, Advanced LIGO, and LISA in Fig. (|4]). The straight solid lines in all plots delineate 
the region of parameter space for which radiation from string loops is in conflict with BBN. In the 
simplest case, this restricts the string tension to satisfy Gji < 7 x 10~^ [63j. Because no loops will 
exist for times after t^, we have shown where this bound can no longer be reliably applied with 
the vertical kink in the solid curve. Interestingly, when loops dissappear before BBN, the tension 
remains constrained by the stochastic radiation from beads. 



where 
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5.3 Stochastic Background Observation and Constraints 
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Figure 5: The upper graph shows the total density of gravitational waves. The value is nearly 
independent of m for a large part of parameter space. We have traced out the portion of the 
parameter space excluded by BBN from loops (solid curve) and beads (dashed curve). Interestingly, 
BBN probes values down to the expected theoretical bound > /i. Because the contours are 
linearly independent of the contours in the burst figures, it may be possible to determine both the 
bead mass and the string tension. On the lower graph we show the spectrum in the LIGO band. 
The sensitivities of Initial and Advanced LIGO are also shown ^64j. 

With some simple approximations, it is possible to determine the main properties of the stochas- 
tic background per logarithmic frequency interval 



^gw(/) 



47r2_ 
3Hl 



nlmax dJi 

dz I dl h'^{f, z,l) 



Imin 



dzdl 



(92) 



The answer does not depend much on the upper bounds of the integrals and they can be taken to 
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Figure 6: Spectrum in the pulsar band. The current bound of 3.9 x 10 ^ is shown [18 



infinity. The lower bounds are set by 
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(93) 



For a decaying network of cosmic strings, all the energy density is transferred to gravitational 
waves. If the decay is relatively sudden, the energy density in the gravitational waves is going to be 
dependent only on the energy in the network before decay, and is independent of the time at which 
the network decays. We therefore expect the stochastic background to be only weakly dependent 
of K, in the case when the network disintegrates during the radiation era. In addition, due to 
the peculiar waveform associated with the ultra-relativistic beads, the spectrum is independent of 
frequency. If the network always obeyed the scaling solution (with energy density pcs oc Gfi/t"^), 
then we would expect the stochastic background to be proportional to Gfi. Because short strings 
scale like matter, there is a deviation from scaling and there is a growth of string energy density for 
times between the short string epoch {t > t^,) and the demise of the network at t^.^, (or tf.q, whichever 
comes first). 

The deviation from scaling between t^, and t** provides an important factor of ^/u^ju = 
_ 1 

(FG/i) 4 to the energy density. Because of this, the stochastic background actually obeys 



(94) 



which is an attractive feature of Fig. ([s]) since it allows much lower tensions to produce an observable 
signal. A similar effect is well known for cosmic string loops[l3]. We computed the total gravitational 
power and spectrum numerically, as shown in Fig. ^ and ([6]). The solid lines in all plots represents 
the region of parameter space for which radiation from string loops is in conflict with BBN. In the 
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simplest case, this restricts the string tension to satisfy Gji < 7 x 10~^ [63]. Again, because no 
loops will exist for times after t^,, we have shown where this bound can no longer be reliably applied 
with the vertical kink in the solid curve. Interestingly, when loops are eliminated before BBN, the 
tension remains constrained by the stochastic radiation from beads. This is illustrated with the 
dashed line in all plots, defined by the BBN constraint VL^^{z > -Zbbn) < 1.5 x 10~^. Here, VL^^ is 
defined as the quantity measured today from all sources more distant than 

'^BBN — 5.5 X 10 . 

Not only are segments of all sizes radiating, but each burst spans a tremendous range of fre- 
quencies, from the length of the segment to well beyond the Planck scale. Because of the scale 
invariance, the total power is simply expressed as 

fi,„ = r],„(/)ln^ (95) 

Jmin 

for any / e [/mm, fmax], where the upper and lower bounds are related to the trans-Planckian and 
Hubble scales respectively, during the short string epoch, as is illustrated in Figs. ([T]). Of course our 
calculation breaks down past Planckian frequencies, and the fact that we have enough energy to 
excite these very high energy gravitons may indicate that massive, universally coupled fields (such 
as the dilaton in string theory) are important. In any case, the ultra- high frequency background is 
totally decoupled from gravitational wave experiments. 




20 40 60 80 20 40 60 



Figure 7: On the left we show the spectrum for = 10~^. The spectrum is scale invariant 
over wavelengths between the Hubble size and well beyond the Planck size. Nearly half of the 
gravitational wave energy density was emitted at super-Planckian frequencies. On the right we 
have the spectrum for Gfi = 10~^^. The IR and UV cutoffs are narrower than for high tension 
strings. From this figure, it is clear that lower frequency experiments, e.g. pulsar timing, are 
ill-suited to probe a wide range of the bead mass, especially for low tension strings. 

It can be understood why the lower tension strings will not produce the Hubble-limited long 
wavelengths found in the higher tension cases (compare the two cases in Fig. [T]). Because the lower 
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string tensions allow a longer period of scaling violation {i.e., segments scaling like matter during 
the radiation era), the Hubble length grows much larger than the typical segment size during their 
radiation peak near t^,^,. Similarly, the lighter strings do not bring the beads to the same speeds 
occurring for higher tension neworks, and the UV cutoff is lower as well. 

6 Conclusion 

The most model independent constraint on cosmic strings, whether meta-stable or stable, comes 
from BBN, which we find restricts the string tension to satisfy (c./. Figjs]) 

G/i < 10-5 (96) 

for any bead mass. This can include strings which decay immediately after formation, and which 
are completely decoupled from further observations. The remaining observables depend upon the 
degree to which cosmic strings are stable. Theories of cosmic strings can only rarely claim the 
strings to be absolutely stable. We find interesting phenomenology for meta-stable strings with 
bead mass within the range 

2 

1 < — < 100 . (97) 

The degree to which this is a fine-tuning depends greatly on the mechanism responsible for string 
breakage. For the bead mass range 40 < m^//i < 80, a stochastic background is detectable by 
Advanced LIGO for tensions > 10^^^. For the longer-lived strings with 70 < iv? / ^ < 80, meta- 
stable strings provide a burst signal detectable by Advanced LIGO for string tensions as low as 
G^x ~ 10~^^. The waveform of these bursts (1//) is very distinctive, and a detection would confirm 
the ultra-relativistic nature of the source. Already, initial LIGO rules out a relatively narrow range 
of K for tensions as low as G^x ~ 10"^". A measurement of these bursts together with the stochastic 
background would enable one to determine both the string tension and the bead mass. This can be 
seen from the linear-independence of the contours in Figs. (|4]&|5]). 

An interesting feature of the predicted stochastic background is the degree to which it is scale 
invariant. For higher tension strings, the spectrum extends between the Hubble scale (at the time of 
emission) to the Planck scale, and possibly beyond. In the upper range, one would expect universally 
coupled massive fields to be radiated with equal power as gravitational radiation [111112]. For lower 
tension strings, the dominant contribution to the gravitational power comes from shorter strings. 
For this reason, low frequency experiments such as pulsar timing are ill-suited to constrain the 
parameter space of meta-stable cosmic strings. 

Eventually, it may also be interesting to seek specific models within string theory where both 
the cosmic string tension and the value of k are predicted and presumably are related (models are 
known where k is much greater than one and is unrelated to the cosmic string tension). We have 
only studied the case of a network of cosmic strings breaking due to beads, although it would be 
very interesting to study the phenomenological signals of the other possible network instability: 
a network of domain walls ending on cosmic strings. We expect that gravitational waves will be 
produced with a similarly characteristic spectrum. 



26 



Acknowledgements 



We would like to thank Jose Blanco-Pillado, Ken Olum, Amir Said, Alex Vilenkin, and Alan 
Wiseman for useful discussions. L.L. would like to thank the KITP where this work was completed. 
This research was supported in part by the National Science Foundation under Grant No. PHY05- 
51164, PHY-0758155, PHY-0505757. The work of B.S. was supported in part by DOE Grant No. 
DE-FG02-91-ER-40672 and NSF Grant No. 0353314. 



A Decay Rate from the Path Integral 

Let us calculate the string decay rate (following [65j). If we treat the string's ground state as an 
approximate energy eigenstate, the decay rate per unit length is equal to twice the imaginary part 
of the energy density. The Euclidean action is given by 

5'_B = /i / dA + m ds 

JY, JdT, 

/•2tt JllfW 1-2-K 



def^^ + m / de^p^{e) + p^{e) 



/o 2 Jo 



where the degree of freedom is the boundary, dT, p{6). We will only consider the portion of S 
that fits in a large rectangle of side lengths X and T, and we assume the string is fiat. The path 
integral gives the vacuum energy density relation 



(0| exp(-WT)lO) = N J [Dp] exp(-5^[p]) (9J 



We define /i by the perturbative result {H) = /x, and the leading corrections to this can be 
calculated in the semiclassical dilute instanton gas approximation. It can be readily shown that the 
above action has saddle points for (9S equal to any number of non-overlaping circles of radius m/p,. 
For a solution with n such holes, the value of the action satisfies 

SE[n X 5^] - fiXT =: nSo = wnm^ / pi . (99) 

In the semi-classical approximation, we expand the path integral around its saddle points. The 
dilute gas approximation is valid in the limit A, T — > cxd, so we ignore the overlap and write the 
path integral equation as 



We are summing over the number of instantons and their positions, dividing by the symmetry 
factor n\. The quantity det' is the determinant of the second functional variation of the action, 
representing fluctuations about a single instanton. (We have set the perturbative contribution to 
the determinant to one by putting the physical parameters into the action.) The prime indicates 
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that wc ignore the zero modes, since these correspond to the position, which has already been 
integrated over. The instanton is expanded in eigenfunctions of the second functional variation as 



cos(i^) i < 



j& ^ y sm{je) j >0 

The functional measure is 



[^^] - n4| = ^-^n^, (102) 



where the prime indicates we remove the zero modes (j = ±1). The single negative mode {j = 0) 
corresponding to dilatations ensures an imaginary contribution to the vacuum energy, i.e., a non- 
zero decay rate. Associated to this negative mode is an overall factor of 1/2, since the path integral 
is divergent only for cq +oo, and not as cq — oo. (This is the factor that Sagredo missed.) 
The result is 



with 



/ -XTe~^° \ 
{exp{-nXT)) = iVexp f -^^-^==- - pXT j (103) 

det' := det' [-/X(9^ - /x] (104) 



We evaluated the determinant using the Riemann zeta function. Hence 



{H) = fJ- — i— exp (— 7rm^//i) , (106) 
and so the bead pair production rate per unit length is 

r2 = ^exp(-7rmV/.) . (107) 

B Calculating the Number Density n{l,t)dl 

The number density of segments satisfies the Boltzmann equation 

^ = -|(MM))-3^n(M)+y, (108) 
where we distinguish between three different source terms: 

9 = 5'loop + 9ic + 5'break (109) 
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with 

gbreakdl = ^2 (^2 ^ n{l' , t)dl' - ln{l ,t)^ dl (110) 

g^^dl = n{l',t)dl' -ln{l,t)^dl (111) 

^loop = (difficult) (112) 

Thus there are three interactions: Loop production is the process by which the network loses 
energy in the form of small loops which fragment and decay. Intercommutation causes two segments 
to trade a quantity of string, keeping the total energy and number of segments fixed. Breakage 
increases the number of segments but leaves the energy unaffected. 

We motivated the breakage formula in the main text. The intercommutation interaction can 
be motivated as follows: A string of length / will intercommute at a rate equal to its length times 
the intercommutation rate per unit length, Tlf . (This rate is easily determined from the scaling 
solution, but will not be important for us.) This explains the second term. The creation of segments 
of length / must involve a string longer than Z/2, since we only need to consider two-string collisions. 
This explains the integral, and the coefficient is over-constrained by the fact that intercommutation 
should leave both the energy and number of segments constant, each of which demands the |. 
Since intercommutation of segments preserves both the total number of segments and their average 
length, we will ignore its effect on n{l,t)dl. (The essential feature of this interaction is it attempts 
to evolve n{l,t)dl into a power-law distribution. Intercommutation will extend the large / tail of 
n{l, t)dl, but this is of little relevance for the short strings which actually radiate.) 

Unfortunately, gioop is difficult to write down, since it depends strongly on how wiggly the strings 
are, which depends on how much string is present. Luckily, we can infer its value from the scaling 
solution. We expect the network will scale as long as most of the string is in superhorizon segments. 
During scaling, two strong effects yield an attractor solution: stretching of superhorizon strings, and 
production of subhorizon loops. Fewer strings cause fewer intercommutations, leading to straighter 
strings and lower loop production. Stretching then causes the string energy density to grow toward 
the attractor value. Too many strings will cause an increase in intercommutation events, making 
strings more wiggly, and more prone to loop production. In the spirit of coarse graining, we will 
think of loop production simply as a continuous shortening of string segments: 

— > /loop (113) 

A key feature is that a segment's length loss due to loop production is proportional to the segment 
length. There are thus three important contributions to /, namely Hubble stretching, gravitational 
radiation, and loop production. 

i = in + /gw + /loop (114) 

Here we will define loop production as being consistent with the scaling solution, whereby the string 
energy density scales like Pcs ~ A*/^^- Because lioop oc /, there is a unique way to do this. 

Gravitational radiation is only important for very short segments. Scaling should cease when 
loop production becomes negligible. This will happen if the strings become smooth, due to reduced 
intercommutation. We expect such behavior during the short string epoch, when the network will 
be dominated by sub horizon sized segments with a density of a few per Hubble volume. Since 
these segments will have slow peculiar velocities, intercommutation will be very rare. 
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Strings longer than twice the Hubble length are always dominated by potential energy, and so 
the length (energy) is increased by expansion: 



;il5) 



For long segments, (f ^) ~ .4, depending on how much intercommutation is going on. If we set 
g = 0, the super- horizon string network would dominate the universe: 



dlfiln{l, t) 



a 



(116) 



We know that the network will obey the scaling solution during the long string epoch. Since a 
scaling network obeys 



Pc 



filn{l,t)dl ~ /i/t^ 



we can use this constraint to find the effective / simply by integrating Eqn.(66). The result is 
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3HI 



t 



(scaling solution) 



(117) 

(118) 
(119) 



As discussed in the text, the scaling solution will apply for t <t^,, after which 

/ ^ -TGfi . (120) 
With an explicit expression for /, we can then solve the Boltzmann equation exactly, obtaining 



n{l,t) 



t<U<t 



for < teg, and 



n{l,t) 



U f2 



(2r2)2e~2r2«t 

(t+t^qf -p2^-r9Mt+Uq)„-(t-U)(t+2teg)/2t'i 

^2 J- 2"^ e 



eq 



^eg — t=¥ t 



;i2i) 



(122) 



for teg < t^,. 



C Gauge Invariant Degrees of Freedom 

In order to obtain the gauge invariant part of the gravitational wave we only need consider the 
spatial components 



„^ I cos 9 — 1/ cos 9 

-T^\k) cos^-l/cos^ 

* cos^ + l/cos^ 








(123) 
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The transverse-traceless part of the polarization tensor can be obtained using the projection 
operator [66] 

p,, = - n.dj, (124) 

as follows 

e^^ = PtP -\Tr{Pe). (125) 
If the wave were to lie entirely on the z-axis the projector would be 




i^j = I 1 I . (126) 
So that some arbitrary polarization tensor 

Cij = I £21 622 £23 I ■ (127) 



en 


ei2 


ei3 


£21 


622 


£23 




£32 


£33 



becomes 

((en -622)72 ei2 

£21 (£22 - £ii)/2 I . (128) 



The so called +- and x -polarisations are defined by 

e+ = |(en-£22) (129) 

and 

£x = £12 = £21- (130) 

For a wave not traveling in the ^-direction we can perform a rotation so that the wave vector 
does lie on the z-axis and then apply the arguments outlined above. This can easily be done noting 
that 

RMRz{\-(t^)^ = ^. (131) 
The rotated metric perturbation polarization tensor is therefore 

e' = RMR^^l - <P)eR^il - <P)RliO), (132) 



with the rotation matrices defined by 



and 



cos a — sin a 

Rzia) =1 sin a cos a | , (133) 
1 



1 

Rx{cf) = I cos a —sin a | (134) 

sin a cos a 
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as usual. 

If we perform the rotation Eq. ( 132 ) on Eq. ( 123 ) we obtain 



/ ^^^^^ 



e^^.(x,^) = —T'^ik) I -^-^ -2sm9 | . (135) 

-2sine We^ 



The transverse traceless part is therefore 



TT / N 2Gm cos^ 9 — 1 ^r.^ , , , i „ _ „ i , , 

ej^ X,.; = ^T°3 0-10 136 

■' cos O 



1 
0-10 




from which the +- and x -polarisations can easily be read off. 

2Gm cos^ e - 1 
r coso' 

and 



T^\k) (137) 



ex = 0. (138) 



D Gravitational Waveform at High Frequency 

Here we compute the gravitational waveform produced by a discontinuity in the bead acceleration 
which can occur when beads cross on a straight string or when a bead encounters a kink on a string. 
We start from the results of Martin and Vilenkin [30] and calculate the gravitational waveform for 
a fly-by of a bead - antibead pair using a Taylor expansion around the time of the crossing. We 



expand xi(t) around t = writing Xi{t) ~ Vot. The Fourier transform of Eq. (10) then reads 

T'^ik) = — / dti^.vo -a\t\) [e'^-' - e^-+*] , 
where uj± = uj{l ± wqcos^). This integral can be done trivially. The integral, 

— / cit(7o^o-a|t|)e-* = 
where we've thrown out a static field piece proportional to 5{uj). This means, 



(140) 



T°3(fc) 



am 



(141) 



Similarly as in the main text, the gravitational wave is linearly polarized and the plus polarization 
can be written as 



2G am 



X,Co' 



27r7r 



1 



cos^ Q - 1 
cos d 



(142) 
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If the beads are relativistic the radiation is beamed in the direction of the bead velocity. If the 
direction of observation is along that direction then one of the l/tu^ or l/cu^ is large and the other is 
0(1), i.e., we are only observing the radiation from one of the two beads. If we choose the direction 
so that cos^ ~ 1, since a = n/mwe can write the waveform as 

2G/i 1 cos^g-1 



In the region where most of the energy is radiated, i.e. Angles 6 ~ I/70 with 70 ^ 1, we can use 
the small angle approximation and set vq ^ 1 to obtain, 

e,(x..)«^m^ (144) 



rv27r7r ^ UJ 

Note that although we have done the computation for the fly-by of a bead - antibead pair the 
waveform is valid up to factors of 0{1) for a bead encounter with a kink on the string which also 
produces a discontinuity in the acceleration. In the case of a fly-by, if the beads miss each other 
with impact parameter b, or for a bead encounter with a kink smoothed on such a length scale, 
we expect the waveform to be suppressed at frequencies above b~^. For a bead - antibead pair 
connected by a string of length / we expect the impact parameter b to be of 0{l), so unless the 
strings are very straight most of the signal will be at frequencies near l^^. On the other hand kinks 
smooth gravitationally very slowly, remaining sharp for the lifetime of the bead - antibead system. 
Therefore beads connected by a kinky string will encounter sharp kinks and radiate every oscillation 
for the lifetime of the system. Neglecting gravitational wave back-reaction, the kink-bead collision 
can be solved exactly. We can write the strain at frequency / for a burst produced at a distance r, 

M/.U)«7r^^f^V (145) 



(1 + Z)(pr{z) \Tnf 

The amplitude is suppressed by an extra factor of {I + z) due to the redshifting of the frequency / 
(see the discussion in section IV of [16]). 
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